Three Dimensional Simulation of Gamma Ray Emission from 
Asymmetric Supernovae and Hypernovae 



O 

o 

(N 



(N 

> 
O 
oo 



in 
o 

o 



Keiichi Maeda 

Department of Earth Science and Astronomy, Graduate School of Arts and Science, University of Tokyo, 

Meguro-ku, Tokyo 153-8902, Japan 



ABSTRACT 

Hard X- and 7-ray spectra and light curves resulting from radioactive decays are computed 
for aspherical (jet-like) and energetic supernova models (representing a prototypical hypernova 
SN 1998bw), using a 3D energy- and time-dependent Monte Carlo scheme. The emission is 
characterized by (1) early emergence of high energy emission, (2) large line-to-continuum ratio, 
and (3) large cut-off energy by photoelectric absorptions in hard X-ray energies. These three 
properties are not sensitively dependent on the observer's direction. On the other hand, fluxes 
and line profiles depend sensitively on the observer's direction, showing larger luminosity and 
larger degree of blueshift for an observer closer to the polar (z) direction. Strategies to derive 
the degree of asphericity and the observer's direction from (future) observations are suggested on 
the basis of these features, and an estimate on detectability of the high energy emission by the 
INTEGRAL and future observatories is presented. Also presented is examination on applicability 
of a gray effective 7-ray opacity for computing the energy deposition rate in the aspherical SN 
ejecta. The 3D detailed computations show that the effective 7-ray opacity ~ 0.025 — 0.027 
cm^ g~^ reproduces the detailed energy-dependent transport for both spherical and aspherical 
(jet-like) geometry. 
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1. INTRODUCTION 

Gamma rays emitted from radioactive isotopes, 
which are explosively synthesized at a supernova 
(SN) explosion (Truran, Arnett, & Cameron 1967; 
Bodansky, Clayton, & Fowler 1968; Woosley, Ar- 
nett, & Clayton 1973), play an important and 
unique role in emission from a supernova, not only 
at 7-ray energies but also at the lower energies: 
from X-rays to even optical or near-infrared (NIR) 
band. Up to a few years after the explosion, the 
decay chain ^^Ni ^ '^^Co ^ ^^Ye (Clayton, Col- 
gate, & Fishman 1969) dominates the high energy 
radiation input to a supernova, with minor contri- 
butions from other radioactivities such as ^''Ni 



^^^Co 5^Fe: Clayton 1974). The decay produces 
7-ray lines with average energy per decay ~ 1.7 
MeV (^^Ni decay with the e-folding time 8.8 days) 
or ^ 3.6 MeV (^^Co decay with the e-folding time 
113.7 days). These line 7-rays are degraded in 
their ways through the SN ejecta, predominantly 
by compton scatterings (see e.g., Casse & Lehoucq 
1990 for a review). Non-thermal electrons pro- 
duced at the scatterings and other processes (pair 
production and photoelectric absorption) heat the 
ejecta, yielding thermal emissions at optical and 
NIR wavelengths. 

The 7-ray emissions from supernovae provide 
a unique tool to study the amount and distri- 
bution of predominant radioactive isotopes ^^Co 
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(therefore ^^Ni produced at the explosion). Den- 
sity structure in the supernova ejecta could also be 
inferred by modeling line profiles that are affected 
by compton scatterings. For example, unexpect- 
edly early detections of the high energy emission 
from SN 1987A (Dotani et al. 1987; Sunyaev et 
al. 1987; Matz et al. 1988) revealed the impor- 
tant role of Ray leigh- Taylor instability and mixing 
in the SN ejecta (e.g., Chevalier 1976; Hachisu et 
al. 1990). The example highlights the very im- 
portance of modeling and analyzing 7-ray emis- 
sion from supernovae. Except the very nearby 
SN 1987A, unfortunately there are to date only 
a few other examples of possible detection of 7- 
rays from the ^^Ni decay chain from supernovae: 
One marginal detection (SN la 1991T: Lichti et 
a. 1994; Morris et al. 1997) and two upper limits 
(SNe la 1986G and 1998bu: Matz & Share 1990; 
Leising et al. 1999). However, now that the In- 
ternational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) has been launched and some new 
gamma ray observatories are being planed (e.g., 
Takahashi, Kamae, & Makishima 2001), it should 
be important to make prediction of 7-ray emis- 
sions (spectra and light curves), taking into ac- 
count recent development of supernova researches, 
i.e., multi-dimensionality (see e.g., Maeda, Maz- 
zali, & Nomoto 2006b and references therein). 

Most of studies on 7-ray transport in SN ejecta 
have been restricted to one-dimensional, spherical 
models. Lately, a few 3D 7-ray transport com- 
putations for SNe have become available: Hoflich 
(2002) developed a 3D 7-ray transport computa- 
tional code and applied it to a 3D hydrodynamic 
model of a Type la supernova (Khoklov 2000). 
Hunger ford, Fryer, & Warren (2003) also devel- 
oped a 3D 7-ray transfer code. They computed 
7-ray spectra for 3D asymmetric (bipolar) type H 
supernova models, and later for 3D single lobe ex- 
plosion models (Hungerford, Fryer, & Rockefeller 
2005). 

Given higher density and smaller ^^Ni mass for 
core-collapse supernovae than SNe la, there is no 
doubt that core-collapse supernovae are more dif- 
ficult to detect in the high energy range (e.g., 
Timmes & Woosley 1997). For example, Hoflich, 
Wheeler, & Khoklov (1998) predicted the detec- 
tion limit of SNe la by INTEGRAL ~ 10 Mpc 
(~ 3 per year), while Hungerford et al. (2003) es- 
timated that for SNe H (more or less similar to SN 



1987A) - 650 kpc (-1-2 per 100 years). Ac- 
cordingly, the number of theoretical prediction of 
hard X-ray and 7-ray emission from core-collapse 
supernovae is still small to date, except models for 
the very nearby SN 1987A (e.g., McCray, Shull, & 
Sutherland 1987; Woosley et al. 1987; Shibazaki 
& Ebisuzaki 1988; Kumagai et al. 1989). Despite 
this, in view of proposed large asymmetry in core- 
collapse supernovae (see e.g., Maeda & Nomoto 
(2003b) and references therein) and its possible 
direct relation to high energy emissions, theoreti- 
cal prediction of high energy emission for a vari- 
ety of core-collapse supernova models is important 
in order to uncover the still-unclarified nature of 
the explosion and understand trends that may be 
possible to observe with current and future instru- 
ments. 

In this respect, potentially interesting targets 
among core collapse supernovae are very energetic 
supernovae, often called "hypernovae" (Iwamoto 
et al. 1998). A prototypical hypernova is SN 
1998bw discovered in association with a gamma 
ray burst GRB980425 (Galama et al. 1998). 
Its broad absorption features in optical spectra 
around maximum brightness and very bright peak 
magnitude led Iwamoto et al. (1998), assuming 
spherical symmetry, to conclude that the kinetic 
energy E^i = i?K/10^^ ergs ^ 30, the ejecta mass 
Mej IOM0, the main sequence mass M^s ~ 
40Mq, and the mass of newly synthesized radioac- 
tivity 5^Ni M(56Ni) ~ O.6M0 (see also Woosley, 
Eastman, & Schmidt 1999). Following this and 
motivated by the deviation between the spheri- 
cal hypernova model prediction and observations 
after r-^ lOO days, Maeda et al. (2006ab) have 
presented comprehensive study comparing various 
observations of SN 1998bw with theoretical expec- 
tations from jet-like aspherical explosion models of 
Maeda et al. (2002), using multi-dimensional radi- 
ation transport calculations. They found that an 
aspherical model with £51 ~ 20 provides a good 
reproduction of optical emission from the explo- 
sion to > 1 year consistently. 

This paper follows the analysis of optical emis- 
sion from aspherical hypernovae applied to SN 
1998bw by Maeda et al. (2006b). In this paper 
we present theoretical predictions of high energy 
emission from the same set of aspherical hyper- 
nova models as presented in Maeda et al. (2002, 
2006ab). Because the models are intrinsically as- 
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pherical, wc make use of fully 3D hard X- and 
7-ray transport computations. In §2, details of 
computational method are presented with a brief 
summary of the input models. In §3, overall syn- 
thetic spectra are presented. §4 focuses on line 
profiles. In §5, light curves of some lines are pre- 
sented. 

In addition to modeling high energy emission 
from SN 1998bw-like hypernovae, also interesting 
is applicability of a gray effective absorptive opac- 
ity for 7-ray transport, which is often used as an 
approximation in computations of optical spectra 
and light curves of supernovae to save computa- 
tional time (e.g., Sutherland & Wheeler 1984). 
Although it has been examined in ID spherically 
symmetric cases and concluded that the approx- 
imation is good if an appropriate value is used 
for the effective 7-ray opacity, it has not been ex- 
amined yet if this is also the case for aspherical 
models. In §6, applicability of the gray absorptive 
7-ray opacity for (jet-like) aspherical models is ex- 
amined by comparing results with detailed 7-ray 
transport and with gray transport. Finally, in §7 
conclusions and discussion, including an estimate 
of detectability of the high energy emission, are 
presented. 

2. METHOD AND MODELS 

2.1. Method 

We have developed a fully 3D, energy-dependent, 

and time-dependent gamma ray transport compu- 
tational code. It has been developed following the 
individual packet method using a Monte Carlo 
scheme as suggested by Lucy (2005). The code 
follows gamma ray transport in SN ejecta dis- 
cretised in 3D Cartesian grids {xi,yj,Zk- linearly 
discretised) and in time steps (t„: logarithmi- 
cally discretised). For the ejecta dynamics, we 
assume homologous expansion, which should be a 
good approximation for SNe la/Ib/Ic. The den- 
sity at time interval (t„, f„+i) is assumed homoge- 
neous in each spatial zone with the value at time 

The transport is solved in the SN rest frame. 
The expansion of the ejecta is taken into account 
as follows. First, cross sections for interactions 
with SN materials given in the comoving frame 
are transformed into the rest frame (Castor 1972). 
The fate of a photon is then determined in the 



rest frame. If a packet survives as high energy 
photons after the interaction, the new direction 
and the energy are given at the comoving frame, 
depending on the specific interaction taking place 
(see below) . The direction and the energy are then 
converted to the SN rest frame (Castor 1972). 
Gamma ray lines from the decay chains ^^Ni 

^ 56Cq ^ 56pg ajj^j 57^1 ^ 57Cq ^ 57pg gj.^ jjj. 

eluded. The numbers of lines included in the com- 
putation are 6 (''^'^'Ni decay), 24 (-"^^Co), 4 ("Ni), 
and 3 (^'^Co) (Lederer & Shirley 1978; Ambwani 
& Sutherland 1988). In the present study we focus 
on the ^^Ni ^^Co ^^Fc chain, which includes 
812 KeV (^'^Ni), 847 keV (^^Co), and 1238 keV 
(56 Co) lines. 

For the interactions of gamma rays with SN 
materials, we consider pair production, compton 
scattering, and photoelectric absorption. For pair 
production cross sections, a simplified fitting for- 
mula to Hubbell (1969) is adopted from Ambwani 
& Sutherland (1988). For compton scattering, the 
Klein-Nishina cross section is used. For photoelec- 
tric absorption, cross sections compiled by Hoflich, 
Khoklov, & Miiller (1992) from Veigele (1973) are 
used. 

If an interaction takes place, the fate of the 
gamma ray packet is chosen randomly in propor- 
tion to the cross section of each possible interac- 
tion. If the pair production is the fate, an elec- 
tron and a positron are created. For the electron 
path, the electron deposits the entire energy to the 
ejecta. For the positron path, the positron anni- 
hilates with an ambient electron, producing two 
7-rays (assuming no positronium formation). The 
positron kinetic energy is deposited to the thermal 
pool. These processes are assumed to take place 
is situ. According to this prescription, a packet 
either becomes the 511 kcV 7-ray lines or is ab- 
sorbed. The fate is selected randomly in propor- 
tion to the branching probability depending on the 
initial photon energy before the pair production. 
If the compton scattering is the fate, the polar an- 
gle of the scattering relative to the incoming pho- 
ton direction in the comoving frame is randomly 
sampled according to the Klein Nishina distribu- 
tion using a standard Monte Carlo rejection tech- 
nique (the Kahn's method). The azimuthal angle 
is randomly selected in the comoving frame. The 
packet is now either a gamma-ray packet or a non- 
thermal electron packet (assumed to deposit the 
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entire energy to the SN ejecta immediately), deter- 
mined randomly according to the energy flowing 
into each branch. If the photoelectric absorption 
is the fate, the entire energy of the packet is ab- 
sorbed. Possible X-ray fluorescence photons are 
not followed in the present simulations, but these 
photons are below the low energy cut-off, making 
negligible effects on the energy range examined in 
this paper. 

For a time advancing scheme, the code has two 
options: The one is the exact time advancing, tak- 
ing into account the time interval at each photon 
flight. The other is no time advancing, assuming 
the flight time through the ejecta is negligible. In 
the present work, we use the exact one for com- 
putations of 7-ray light curves (and for computa- 
tions of optical emission in Maeda et al. 2006b), 
while use the no-time advancing for computations 
of spectra. This is in order to optimize the num- 
ber of photons escaping at a given time interval in 
the spectral computations, since the computation 
of spectra with fine energy bins already needs a 
large number of packets and large computational 
time. For 7-rays, the negligible time delay approx- 
imation is a good one. We have confirmed this for 
the present models, by performing the fully time- 
dependent computation but with the number of 
packets entering each time interval smaller than 
used in the standard no-time advancing spectrum 
computations. 

In view of the recent investigation by Milne et 
al. (2004) that not all the published ID gamma 
ray transport codes give mutually consistent re- 
sults, we test the capability of our new code 
by computing gamma ray spectra based on the 
(spherical) SN la model W7 (Nomoto, Thiele- 
mann, & Yokoi 1984), for which many previous 
studies arc available for comparison. Wc have 
compared our synthetic gamma ray spectra at 25 
and 50 days with Figure 5 of Milne et al. (2004). 
We found an excellent agreement between our re- 
sults and the spectra resulting from majority of 
previous codes, e.g., of Hungerford et al. (2003). 

In this study, the ejecta arc mapped onto 60^ 
Cartesian grids. For spectrum synthesis, 1.5 x 10^ 
photon packets with equal initial energy content 
are used. In escaping the eject, these packets are 
binned into 10 angular zones with equal solid angle 
from ^ = 0° to 180° (here is the polar angle from 
the ^;-axis) and into 3000 energy bins up to 3 MeV 



with equal energy interval 1 kcV. For 7-ray light 
curve computations, 10^ photon packets are used. 
In escaping the ejecta, the packets are binned into 
36 time step logarithmically spaced from day 5 to 
day 300, as well as into the angular and energy 
bins. 

2.2. Models 

Input models for the 7-ray transport computa- 
tions are taken from the aspherical model A and 
the spherical model F from Maeda et al. (2006b). 
Model A is a result of a jet-like explosion with the 
initial energy input at the collapsing core injected 
more in the z-axis than in the r-axis (see Maeda 
et al. 2002 for details). In Maeda et al. (2006ab), 
we examined optical light curves (< 500 days) as 
well as expected optical spectral characteristics in 
both early (< 100 days) and late (> 100 days) 
phases. The structure of Model A at homologous 
expansion phases is shown in Figure 1. 

In the present work wc examine the follow- 
ing models: {M^-JMq, E51, M(56j^i)/M0) = 
(10.4,10,0.31) and (10.4,20,0.39) for Model A, 
and (10.4,10,0.28) and (10.4,50,0.40) for Model 
F (hereafter Mej is the ejecta mass, E51 is the 
kinetic energy of the expansion in 10''^ ergs, and 
M(^^Ni) is the mass of ^^Ni synthesized at the 
explosion). Maeda et al. (2006b) concluded that 
optical properties of SN 1998bw arc explained con- 
sistently by Model A with the energy E51 ~ 20 
and with M(^^Ni) ~ O.4M0, so that we regard 
Model A representing a prototypical hypernova. 

As seen in Figure 1, Model A is characterized 
by concentration of ^^Ni distribution along the z- 
axis, which is a consequence of the explosive nu- 
cleosynthesis in jet-like aspherical explosions (Na- 
gataki 2000; Maeda et al. 2002; Maeda & Nomoto 
2003b). Because ^^Ni is the source of 7-rays, the 
distribution will affect the 7-ray transport and re- 
sulting hard X- and 7-ray emissions sensitively. 
Figure 2 shows the distribution of ^^Ni along the 
line of sight for our models. The amount of ^^Ni is 
integrated in the plane perpendicular to the line of 
sight within the constant line-of-sight velocity in- 
terval. Such that, Figure 2 shows a profile of 7-ray 
lines from the decay of ^^Ni or ^^Co in optically 
thin limit. 

Figure 2 shows that Model A yields the '"'^Ni at 
high velocities if it is viewed from the polar (z) 
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direction. On the other hand, if it is viewed from 
the equatorial (r) direction, it yields the distri- 
bution sharply peaked at zero and low velocities. 
Note that Model A with £^51 = 20 shows ^^Ni at 
velocities higher than in Model F with E^i = 50. 

3. HARD X- AND 7-RAY SPECTRA 

Figures 3 and 4 show synthetic 7- and hard X- 
ray spectra at 25 days after the explosion. Figures 
5 and 6 show the ones at 50 days after the explo- 
sion. Thanks to the largo Sx/A^ej as compared to 
(normal) SNe II and to the absence of a massive 
hydrogen layer, the emergence of the high energy 
emission is much earlier for the present models (in 
the order of a month) than SNe II (in the order of 
a year). 

An aspherical model yields the emergence of 
the high energy emission earlier than a spherical 
model with the same energy: At 25 days, Model A 
with E51 = 10 already shows high energy appear- 
ance although Model F with Er^i = 10 docs not 
show up. Indeed, the former, despite the energy 
£^51 = 10, gives the date of emergence comparable 
to Model F with £51 = 50. The time scale of the 
7-ray emission will be further discussed in §5. 

Model A is characterized by the large line-to- 
continuum ratio, especially at early on (Figure 
3). Because the continuum is formed by degra- 
dation of line photons by compton scatterings, 
the large line-to-continuum ratio is realized for 
the plasma with low optical depth for compton 
scatterings. Indeed, 7-ray transport computations 
typically predict larger the ratio at more advanced 
epochs (see e.g., Sumiyenov et al. 1990). It is also 
seen by comparing Model F with different energies 
(£51 = 10, 50) at day 50 (Figure 5): The larger en- 
ergy leads to smaller optical depth, and therefore 
to the larger line-to-continuum ratio. At 25 days, 
Model A with either £51 = 20 or 10 yields the ra- 
tio larger than Model F with £51 = 50, which is 
attributed to the existence of high velocity ^^Ni at 
low optical depth in Model A (Figures 1 & 2). At 
50 days, the effect becomes less significant (i.e., the 
ratio becomes comparable for Model A (£51 = 10 
or 20) and for Model F with £51 = 50), but still 
visible as compared with Model F with £51 = 10 
(Figure 5). 

Another feature is seen in hard X-ray spectra. 
Model A has the hard X-ray cut-off at energy 



higher than Model F (Figures 4 and 6). The cut- 
off is formed by photoelectric absorption, which 
is dependent on metal content (e.g., Grebenev & 
Sunyacv 1987). At 25 days, only the emission 
near the surface (along the 2-axis for Model A) 
can escape out of the ejecta. The cut-off is there- 
fore determined by metal content near the surface. 
In Model F, the surface layer is dominated by in- 
termediate mass elements (i.e., a CO core of the 
progenitor star). On the other hand, in Model 
A the emitting region contains a large fraction of 
Fe-peak elements, most noticeably ^^Ni (or Co, 
Fc). Therefore, the photoelectric cut-off should 
be at the energy higher in Model A than Model F. 
At more advanced epochs, an observer looks into 
deeper regions. Then the contribution from the 
deeper region becomes bigger and bigger, yield- 
ing increase of the cut-off energy in Model F since 
'''^Ni is centrally concentrated (Figure 2). This ef- 
fect is also seen in Model A, but to the smaller 
extent. Although the mass of ^^Ni within a given 
velocity interval increases toward the center (ex- 
cept the inner most region where the ^^Ni fraction 
is very small) also in Model A along the ^-axis, the 
increase is less dramatic than model F (e.g., com- 
pare the masses of ^^Ni at 20,000 and 10,000 km 
s~^ for Model A in Figure 2). Therefore, the metal 
content of the emitting region does not temporally 
evolve very much in Model A. 

In the above discussion on the photoelectric ab- 
sorption, one would expect to use the ^^Ni distri- 
bution along the "r-axis", not the 2:-axis, for the 
observer at the r-axis. However, it is not the case. 
As discussed by Hungerford et al. (2002, 2005), 
even for the observer at the r-axis, ^^Ni contribut- 
ing the most of the emission is that at the 2;-axis, 
i.e.. an observer at the r-axis looks at the emit- 
ting polar ^^Ni blobs sideways (see also Maeda et 
al. 2006b). To clarify this, in Figure 7 we show 
last scattering points of hard X- and 7-rays. Also 
shown are contours of the optical depth for ob- 
servers at +z- and -|-r-directions. The observer 
at the z-axis looks at the emitting blob moving 
toward the observer, while the observer on the r- 
plane looks at a pair of the emitting blobs side- 
ways. 

This is further discussed in §4 and 5, but here 
we point out one difference related to effects of the 
viewing angle in our model from Hungerford et al. 
(2002). Their models do not yield large difference 
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in the absolute flux, while our models do show 
boost of high energy luminosity toward the z-axis. 
A similar behavior is also seen in optical emission 
in our models (Maeda ct al. 2006b). The mecha- 
nism of the boost should be the same as the one 
for optical photons. Initially at high density only 
the sector-shaped region along the z-axis can yield 
escaping photons (Fig. 7). The cross sectional 
area of this photosphere is larger toward the z- 
axis than the r-axis, making the luminosity boost 
toward the z-axis. As time goes by, the ejecta opti- 
cal depth decreases, and therefore the photosphere 
moves to cover the equatorial region. The differ- 
ence becomes less and less significant, as seen by 
comparing Figures 3 and 5. The different behav- 
ior is probably due to the different degree of the 
penetration of ^^Ni into the outer layers. Hunger- 
ford et al. (2002, 2005) presented SN II models 
with a massive hydrogen envelope, which does not 
exist in our model SNe Ib/c. The existence of the 
hydrogen envelope generally yields less aspherical 
distribution of ^^Ni and the emitting region in SNe 
II than SNe Ib/c (e.g., Wang et al. 2001). In any 
Hungerford et al. (2005) suggested, the 
high energy emission can very sensitively depend- 
ing on the type and degree of asymmetry, therefore 
transport calculations are very important. 

4. LINE PROFILES 

We now turn our attention to line profiles. Fig- 
ure 8 shows line profiles of the ^^Ni 812 keV and 
^''Co 847 keV lines at 25 days. Figure 9 shows the 
5^ Co 1238 keV line. Figures 10 and 11 show the 
same lines at 50 days. At these early epochs, the 
ejecta are not thin to high energy photons (Fig- 
ure 7), and therefore the line profiles are different 
from the ones expected in optically thin limit (Fig- 
ure 2). 

At day 25, an observer only sees the region near 
the surface (Figure 7). For Model A, the region 
is further concentrated along the 2-axis (the top 
of the "^^Ni bubble" in Figure 1). Figures 8 and 
9 show that the line profiles are very asymmet- 
ric, showing only the emission at the blue com- 
ing from the region moving toward an observer, 
except an observer at the r-axis. For an ob- 
server at the r-axis, even the emission at the rest 
wavelength is seen, which is especially evident for 
the more energetic, therefore less optically thick. 



model (Model A with E^i = 20). The emitting 
region is now at the outer edge(s) of the ^^Ni dis- 
tribution, and the density is relatively small there 
(Figure 1). The line-of-sight, connecting the emit- 
ting region(s) and an observer at the r-direction 
passes only through the low density region, includ- 
ing the region moving away from the observer. On 
the other hand, the line-of-sight for an observer at 
the 2;-axis passes through the dense, central region. 
For example, this behavior is seen in Figure 7 by 
comparing the regions, having the optical depth 

less than 10 at day 25, for the observers at 
z- and r-axes. For the observer at +z direction, 
the region contains only the emitting blob moving 
toward the observer. For the observer at the +r 
direction, on the other hand, a pair of the emitting 
blobs along the z-axis, seen from the side by the 
observer, are included in the regions with Tj < 10. 

At 50 days, the line becomes broad redward as 
a photon even from the far side escape out of the 
ejecta more easily because of decreasing density. 
Still, only the blue part is seen for an observer at 
the 2;-direction (Figures 10 & 11). It is in contrast 
to Model F with E51 = 50, which now show the 
emission at the rest wave length. It is interest- 
ing to see that the line profiles of Model A (either 
E51 = 10 or 20) viewed at the ^-direction resemble 
to the ^^Ni distribution in the hemisphere mov- 
ing toward the observer (Figure 2). This suggests 
that the photons from the far hemisphere are al- 
most totally blocked by the high density central 
region, while the ^^Ni-rich region itself is nearly 
optically thin. It is indeed the case as seen in 
Figure 7. For Model A viewed at the r-direction 
and for Model F, the line profile is explained by 
continuously decreasing escape probability. These 
arguments imply that, assuming that wc have a 
temporal series of 7-ray observations, the line pro- 
files for Model A viewed at the 2-direction should 
show the time interval within which the line pro- 
files are almost fixed, while Model A viewed at 
the r-direction and Model F should show the con- 
tinuous changes in the line profile until the entire 
ejecta become optically thin. This could be an 
interesting observational target, since the tempo- 
rally "fixed" line profile, if observed, suggests that 
the viewing angle is close to the pole. In this ease, 
the line shape directly traces the distribution of 
^^Ni. In addition, line profiles, sensitively depen- 
dent on the viewing angle as well as the degree of 
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asphericity, could be used as a tracer of the distri- 
bution of ^^Ni and density once 7-ray observations 
at ~ 1 MeV with sufficient sensitivity become pos- 
sible. The detectability will be discussed in §7. 

5. LIGHT CURVES 

Figures 12 shows light curves of the ^^Ni 812, 
^^Co 847, and '^^Co 1238 kcV lines at the distance 
10 Mpc. For the 1238 keV line, the flux is com- 
puted by integrating a spectrum in 1150 - 1340 
keV range. For the 812 and 847 keV hues, the 
sum of the fluxes is computed first by integrating 
a spectrum in the 800 - 920 keV range, then the in- 
dividual contribution is computed with the decay 
probability at the corresponding epochs assuming 
that the escape fractions for these two lines are 
equal. The assumption of the equal escape frac- 
tions should be a good approximation, since these 
two lines are close in the energy, yielding almost 
identical compton cross sections (see also Milne et 
al. 2004). 

Model A shows the emergence of the 7-rays 
earlier, even though the energy is smaller, than 
Model F. This is due to the large amount of ^^Ni 
at low density and high velocity regions. Later 
on, the 7-rays become stronger for Model F with 
^51 = 50 than for Model A {E51 = 10, 20). Three 
effects are responsible: (1) The contribution of the 
dense central region stopping the photons becomes 
large. (2) The overall optical depth r cx M'^/E is 
larger for Model A than Model F. (3) Model F 
with E51 = 50 has a bit large M(^^Ni) as com- 
pared to Model A. Because of this temporal be- 
havior, the enhancement of 7-ray flux is more no- 
ticeable for the ^^Ni line (with the e- folding time 
8.8 days) than the ^^Co lines (with the e-folding 
time 113.7days). 

As discussed in §3, the effect of the viewing an- 
gle is seen at early epochs. Model A emits more 
toward the z-direction than the r-direction. This 
behavior is qualitatively similar to that in opti- 
cal emission (Macda ct al. 2006b). The effect 
virtually vanishes around ~ 100 days for Model 
A with E51 = 20, which is consistent with a 
rough estimate that the cjccta become optically 
thin to compton scatterings at ~ ^yTQ ~ (1250 x 
iMcj/MQf/E5iy/'^ ~ 80 days (See e.g., Maeda 
et al. 2003a: sec also Figure 7), where tq is the 
optical depth to compton scatterings at day 1. In 



the above estimate, the cross section is assumed 
to be 1/6 of the Thomson scattering cross section, 
and composition Yg = 0.5. 

6. COMPARISON BETWEEN DETAILD 
AND GRAY TRASNPORT 

In previous sections, we presented model pre- 
dictions for high energy spectra and light curves. 
However, 7-ray transport in the SN ejecta has an- 
other role which is as important as the high en- 
ergy emission itself: the 7-rays give heating of 
the ejecta, i.e., non-thermal electrons scattered 
off from ions by compton scatterings rapidly pass 
their energies to thermal particles through ioniza- 
tion, excitation, and scatterings with thermal elec- 
trons. The thermal energy is then converted to 
optical photons. In this way, the energy lost from 
the high energy photons determines the optical lu- 
minosity from a supernova. In this section, we 
examine how much energy is stored in the ejecta. 

In Maeda et al. (2006b), we made use of the 
detailed 3D high energy photon transport to com- 
pute optical light curves and nebular spectra of 
Models A and F (and other models). On the 
other hand, a gray effective absorptive assump- 
tion for 7-ray transport has been frequently used 
for computations of 7-ray deposition and optical 
light curves. Sometimes it has been used, although 
in spherically symmetric models, to compute opti- 
cal light curves for (hypothetical) non-spherically 
symmetric supernovae (e.g., Tominaga et al. 2005; 
Folatelli et al. 2006). It is therefore important 
to examine applicability of the assumption for as- 
pherical supernovae. 

Figure 13 shows synthetic optical light curves 
for Models A and F. Results with the detailed 3D 
transport and with the simplified gray absorptive 
7-ray opacity (with various effective opacity Kj) 
arc compared. Figure 13 shows that the gray ab- 
sorptive approximation is actually good for com- 
putations of 7-ray deposition, thus for computa- 
tions of optical light curves, if an appropriate value 
is used for the effective opacity. With the value 

= 0.027 cm^ g~^, we obtain optical light curves 
for both the spherical model F and the asphcrical 
model A almost identical to those obtained with 
the detailed 7-ray transport computations. 

At late epochs > 100 — 200 days (depending 
on models) the ejecta become optically thin to 7- 
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rays. At these epochs, 7-rays suffer at most only 
one compton scattering before escaping the ejecta. 
In this idealized situation, the "effective" 7-ray 
opacity can be computed by taking an appropri- 
ate average of the Klein-Nishina cross section for 
various scattering angles and 7-ray line energies. 
In this way, Sutherland & Whcclcr (1984) yielded 
K-y = 0.022 cm^ g~^ assuming the typical line en- 
ergy ~ 2 MeV. This value is dependent on the line 
list. Wc have also performed the estimate of the 
effective opacity in the optically thin limit under 
the condition that the averaged absorbed energy 
per photon flight is equal for the detailed case and 
for the effective absorption case, and found that 
for the 56Co lines, = 5.08 (Z/A) cm^ g-i. This 
yields ^ 0.025 cm^ g~^ for the SNc Ib/c com- 
position. Indeed, Figure 13 implies that this value 
is probably even better than 0.027 cm^ g~^ in the 
late phase, although the difference is small. We 
emphasize that in this situation, the effective opac- 
ity is independent from geometry of the ejecta. It 
is confirmed, although not generally, by the fact 
that the n-y = 0.027 cm^ g~^ reproduces the light 
curve obtained by the detailed computations for 
both Models A and F. 

Somewhat surprising is that (1) the same value 
= 0.027 cm^ g~^ reproduces the detailed trans- 
port computations rather well at early epochs, and 
(2) this applies not only to the spherical model, 
but also to the aspherical model. Because mul- 
tiple scatterings take place, now the effective 7- 
ray opacity can be geometry- and time-dependent, 
and can be different from the value at the late 
phases. For example, Sutherland & Whcclcr sug- 
gested K-y = 0.03 cm^ g~^ taking into account mul- 
tiple scatterings. Frannson (1990) gave the value 
K-y — 0.03 cm^ g~^ with YJ, = 0.5 from a series 
of Monte Carlo simulations. Colgate, Petschek, & 
Kriese (1980) gave the value K-y = 0.028 cm^ g~^. 
These authors used different ejecta models. The 
small differences among these studies (although 
there are some exceptions, e.g., 0.07 cm^ g~^ by 
Woosley, Taam, & Weaver (1986)) suggest that 
the effects of density and ^^Ni distribution on the 
total deposition rate is small (in spherically sym- 
metric case) . The present study suggests that this 
is also the case even in asymmetric cases. 



7. CONCLUSIONS AND DISCUSSION 

In this paper, we first reported the develop- 
ment of a new code for 3D computations of hard 
X-ray and 7-ray emission resulting from radioac- 
tive decays in supernovae. The code is useful for 
many problems. Among them is applicability of a 
gray absorptive approximation in supernova ejecta 
without spherically symmetry. We have shown, 
although only for some specific models, this is ac- 
tually a good approximation. This approximation 
has been extensively used even in analyzing super- 
novae with hypothetical asymmetry (e.g., Tomi- 
naga et al. 2005; Folatelli et al. 2006), and we 
have confirmed the applicability using the 3D com- 
putations. 

With the code, we have presented predictions of 
hard X-ray and 7-ray emission in aspherical hyper- 
nova models. The models were verified by fitting 
optical properties of the prototypical hypernova 
SN 1998bw (Maeda et al. 2006ab), therefore we 
regard the model prediction being based on "real- 
istic" hypernovae. 

Irrespective of the observer's direction, the as- 
pherical models yield (1) early emergence of the 
high energy emission, and therefore the large peak 
flux, (2) large line-to-continuum ratio, and (3) 
large cut-off energy in the hard X-ray band, as 
compared to spherical models. These are qualita- 
tively similar to what are expected from extensive 
mixing of '"'^Ni (see e.g., Cassc & Lehoucq 1990) 
by e.g., Ray leigh- Taylor instability (e.g., Hachisu 
et al. 1990). However, degree of the ^^Ni penetra- 
tion is more extreme, and therefore the effects of 
^^Ni mixing is more noticeable, in our models than 
spherical mixing cases. For example, for Model A, 
the cut-off energy (by photoelectric absorptions) 
is already at ~ 100 keV at the emergence of the 7- 
rays, and it does not temporally evolve very much, 
because the mass fraction of ^^Ni (Co) does not 
increase very much toward the center along the z- 
axis. This constant cut-off energy would be useful 
to distinguish the aspherical models from spheri- 
cal ones, once we have a sequence of hard X-ray 
observations. 

Another interesting result is that line profiles 
are sensitively dependent on asphericity, as well 
as observer's direction. Lines should show larger 
amount of blueshift for an observer closer to the 
polar {z) direction. For an observer at equatorial 



8 



(r) direction, lines show emission at the rest wave- 
length even at very early epochs. These are not ex- 
pected in spherically symmetric models. Although 
analysis of line profiles will need future observato- 
ries with very high sensitivity (e.g., Takahashi et 
al. 2001) and/or a supernova very nearby, once 
it becomes possible then it will be extremely use- 
ful to trace asphericity and energy, therefore the 
explosion mechanism, of supernovae. 

For example, a possibly interesting observa- 
tional strategy is implied. In Model A the ^^Ni re- 
gion becomes optically thin at a rather early phase 
when the inner regions are still optic;ally tliic;k. It 
will lead to the following temporal evolution for an 
observer at the ^;-direction. First, a line becomes 
redder and redder with time. Then at some epoch 
it stops changing the shape very much, present- 
ing closely the ^^Ni distribution in the hemisphere 
moving toward the observer. Finally, the line from 
the other hemisphere appears. This evolution is 
unique as compared with spherical models or for 
an observer at r-direction, either of which should 
show continuous reddening until entire ejecta be- 
come optically thin. 

Because we show that there are features of the 
aspherical models, some of which depend on the 
degree of asphericity but not on the observer's di- 
rection, while the others depend on both the as- 
phericity and the direction, combination of vari- 
ous analyses will be helpful to distinguish models 
and observer's directions. A problem is, of course, 
how many hypernovae are expected to be able to 
be observed in current and future observatories. 
For the aspherical model A with E^i ~ 20, the 
expected maximum line fluxes are about a half of 
those of typical SN la models. The less energetic 
model (i?5i = 10) yields even smaller fluxes. This 
is because of smaller amoimt of ■''^Ni and larger ra- 
tio Mej^/i?K (and therefore later emergence of the 
radioactive decay lines) than the SN la models. 
Even worse, the occurring rate is much smaller for 
hypernovae than SNe la. In this respect, there is 
no doubt that SNe la are most promising targets 
in radioactive 7-rays. 

However, hypernovae could still be an interest- 
ing targets in 7-rays among core-collapse super- 
novae. Taking SN 1987A as a typical SN II, its 
peak 847 keV line flux was ~ 10"'^ photons cm~^ 
s~^, yielding 3 x 10~^ cm~^ s~^ if it were at 
10 Mpc. Our aspherical hypernova models predict 



the peak flux more than two orders of magnitude 
larger than this. Therefore, for a given sensitiv- 
ity of instruments, it is visible at least an order 
of magnitude farther than a typical SN II. Low 
mass SNe Ib/c should also be discussed. Although 
the nature of the SNe Ib/c seems very heteroge- 
neous, let us assume that Mcj = 2Mq, £'51 = 1, 
and M(5^Ni) = O.OTM©. These values give the 7- 
ray escaping time scale comparable to the present 
models (see §5), resulting in the estimate of the 
peak line flux ~ 5 times smaller than the present 
models. For the sensitivity of the INTEGRAL, 
we estimate that the maximum distances within 
which 7-rays are detectable are ~ 500 kpc, 4 Mpc, 
and 7 Mpc, for SNe II, Ib/c, and hypernovae, re- 
spectively (here we assume the line width 10, 20, 
and 40 keV for SNe II, Ib/c, and hypernovae, re- 
spectively). The latter two, especially the detec- 
tion limit for hypernovae, cover some star burst 
galaxies like M82 (~ 4Mpc). Since SN 1987A, the 
nearest SNe for each type are SN II 2004dj (~ 3 
Mpc), SN Ic 19941 7 Mpc), and a hypernova 
(but weaker than SN 1998bw) SN Ic 2002ap (~ 10 
Mpc). the prototypical hypernova SN 1998bw was 
at ^ 36 Mpc. Therefore, the sensitivity of the 
INTEGRAL is unfortunately not enough to de- 
tect these core-collapse supernovae, and we will 
have to wait for next generation hard X- and 7- 
ray telescopes. If these telescopes will be designed 
to archive sensitivity better than the INTEGRAL 
by two orders of magnitudes (see e.g., Takahashi et 
al. 2001), then the maximum distances are 5 Mpc, 
40 Mpc, and 70 Mpc for SNe II, Ibc, and hyper- 
novae, respectively. The detection limits for SNe 
Ib/c and hypernovae then cover some clusters of 
galaxies like Virgo (~ 18Mpc), Fornax (~ ISMpc), 
and even Hydra (~ 40Mpc). This sensitivity will 
lead to comprehensive study of 7-ray emission, and 
therefore of explosive nucleosynthesis, in core col- 
lapse supernovae. 

Observed occurring rate is ~ 7 x 10"'^ and 
~ 1 X 10~^ per year in an average galaxy for SNe 
II and SN Ib/c, respectively. The hypernova rate 
is rather uncertain. Podsiadlowski et al. (2004) 
gave a conservative estimate of the rate ~ 10~^ per 
year. These values give, by multiplying the rate 
and the volume detectable, the observed probabil- 
ity ~ 1 (SNe II) : 70 (SNe Ib/c, but not hyper- 
novae) : 4 (hypernovae) for a given detector. The 
difference between usual SNe Ib/c and hypernovae 
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may (probably) be even smaller, since the estimate 
of SNe Ib/c flux above is probably the upper limit 
(i.e., Mcj = 2Mq seems to be almost the lower 
limit) and that the hypernova rate is possibly the 
underestimate. 

In sum, as long as the detectability is con- 
cerned, a hypernova is worse than (usual) SNe 
Ib/c, but potentially better than SNe II. Another 
question is if any asphericity similar to that de- 
rived for the prototypical hypernova SN 1998bw 
(Maeda et al. 2006ab) exists in (usual) SNe Ib/c. 
If it does, as suggested by e.g., Wang et al. (2001), 
then most of the properties shown in this paper 
should also apply for hard X- and 7-ray emission 
from these SNe Ib/c (since the time scale is simi- 
lar to hypernovae: sec §5). One should of course 
take into account the fact that the expansion ve- 
locity and M(^^Ni) are smaller than hypernovae, 
and ultimately needs direct computations based 
on a variety of explosion models. 
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Fig. 1. — Model A (Maeda et al. 2002) with E51 = 20. The left panel shows density distribution in a 
logarithmic scale at 10 days after the explosion. The ejecta is already in a homologous expansion phase, so 
that the distribution is shown in the velocity space. The right panel shows mass fractions of ^^Ni in a linear 
scale. 
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Fig. 2.— Line-of-sight distribution for Model F with E51 = 10 (gray), F with E51 = 50 (black), and 
for Model A with E^i ~ 20 along the polar (z) direction (red) and along the equator (r) direction (blue). 
Masses of ^^Ni within the line-of-sight velocity V ^ V -\- dV (with dV constant) are shown as a function of 
the line-of-sight velocity. The mass is in an arbitrary unit. 
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Fig. 3. — Synthetic hard X-ray and 7-ray spectra at 25 days after the explosion at the reference distance 10 
Mpc. The spectra are shown for Model A ((a,b) E51 = 10; (c,d) E^^i — 20). The observer's direction is along 
the z-axis (a, c) and along the r-axis (b, d). In each panel, a synthetic spectrum for Model F {E51 = 50) 
at 25 days is also shown (gray) for comparison. Model F with E^i = 10 is still extremely optically thick in 
these wavelengths at this epoch, and therefore not visible. 
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Fig. 4. — Same as Figure 3 but with the horizontal axis in logarithmic scale, for presentation of the hard 
X-ray spectra. 
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Fig. 5. — Synthetic hard X-ray and 7-ray spectra at 50 days after the explosion at the reference distance 
10 Mpc. See the caption of Figure 3. Here, two spherical models (Model F with E51 = 10 (lower) and 50 
(upper)) are shown for comparison. 
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Fig. 6. — Same as Figure 5 but with the horizontal axis in logarithmic scale, for presentation of the hard 
X-ray spectra. 
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Fig. 7. — Last scattering points of 7 and hard X-ray photons at 25, 50, 100, and 200 days since the explosion 
for Model A with E51 = 20. Also shown are contours of optical depth for an observer at ^-axis (at +z 
direction: solid) and at r-axis (at +r direction: dotted). The horizontal and vertical axes are respectively 
the Vx and the Vz axes. At 25 days, the contour is shown for the optical depth = 1, 10, and 30, while for 
other epochs it is shown only for = 1. The optical depth for the contour is here computed assuming the 
opacity equal to 1/6 of the Thomson cross section. 



18 



u 1E-4 



S 

'a, 1E-5 



u 

(A 

c 
o 
4^ 
o 
f 
Q. 



1E-6I — ' 
0.80 



(a) Model A (z), E^^=10 
@ 25 day 




0.84 0.88 
Energy (MeV) 



0.92 



'> 1E-4 



0.80 



(c) Model A (z), E,^=20 
@ 25 day 




0.84 0.88 
Energy (MeV) 



0.92 



^ 1E-5 



(A 
C 

o 

o 1E-6 
a. 





1 1 1 1 1 1 1 1 
(b)ModelA(r), E =10 ; 






@ 25 day 









0.80 



0.84 0.88 
Energy (MeV) 



0.92 



o 1E-4 

Q. 

s 



1> 

0) 



1E-5 



O 

o 
f 
Q. 



1E-6I — ' 
0.80 



(d) Model A (r), E^,=20 
@ 25 day 




0.84 0.88 
Energy (MeV) 



0.92 




Fig. 8. — Line emission of ^^Ni 812 KeV and ^°Co 847 KeV at 25 days after the explosion at the reference 
distance 10 Mpc. Models spectra are shown for Model A with E51 = 10 (a,b) and £^51 = 20 (c,d), with the 
observer at the 2;-axis (a,c) and at the r-axis (b,d). Also shown is a spectrum of Model F with £^51 = 50 (e). 
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Fig. 9. — Line emission of ^®Co 1238 KeV at 25 days after the explosion at the reference distance 10 Mpc. 
See the caption of Figure 8. 
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Fig. 10. — Line emission of ^^Ni 812 KeV and ^^Co 847 KeV at 50 days after the explosion at the reference 
distance 10 Mpc. Models spectra are shown for Model A with E51 = 10 (a,b) and -E51 = 20 (c,d), with the 
observer at the 2;-axis (a,c) and at the r-axis (b,d). Also shown is a spectrum of Model F with £^51 = 50 (e) 
and £51 = 10 (f). 



21 




(c) Model A (z), £^,=20 
@ 50 day 



1E-4 



0) 



E 
u 

(A 

o 
o 

Q. 



1E-5 




1.16 1.20 1.24 1.28 
Energy (MeV) 



1.32 



u 
a. 



(e) Model F, E^^=50 
@ 50 day 




1.16 1.20 1.24 1.28 
Energy (MeV) 



1.32 



1E-4 



0) 



E 
o 

(A 
C 

o 
o 

a. 



1E-5 



(d) Model A (r), E^,=20 
@ 50 day 




1.16 1.20 1.24 1.28 
Energy (MeV) 



1.32 



o 1E-5 
cl 



1E-6 



E 
o 
in 
c 
o 

o 
.c 
a 



1E-7 



(f) Model F, =10 

@ 50 day 




1.16 1.20 1.24 1.28 
Energy (MeV) 



1.32 



Fig. 11. — Line emission of ^^Co 1238 KeV at 50 days after the explosion at the reference distance 10 Mpc. 
See the caption of Figure 10. 
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Fig. 12.— Light curves of (a) the ^^^Ni 812 KeV line, (b) the ^^Co 847 KeV hne, and (c) the ^^Co 1238 KeV 

Hnc at the reference distance 10 Mpc. For the 812 & 847 KeV hnes, the flux is computed by integrating a 
spectrum in the 800 — 920 KeV energy range, then by assuming that the escape fractions for these two lines 
are equal. For the 1238 line, the flux is computed by integrating a spectrum in 1150 — 1340 KeV energy 
range. Shown here are Model A with E^i ~ 20 (thick-black) and 10 (thin-black), and Model F with E^i = 50 
(thick-gray) and 10 (thin-gray). For Model A, the flux is shown for the z-direction (solid), the r-direction 
(dotted), as well as the angle-averaged one (dashed). 




Fig. 13. — Optical light curves computed for (a) Model F with £^51 = 50 and (b) Model A with £^51 = 20 
(averaged over all the solid angles). Optical photon transport is solved with a simplified gray representation 
for optical opacity (Chugai 2000) by the Monte Carlo radiation transfer method (Maeda et al. 2006b) . For 
the gamma ray transport, the detailed transport (gray) is compared with the simplified gray atmosphere 
approximation with the effective absorptive gamma ray opacity = 0.13, 0.27, and 0.54 cm^ g^^ (black 
curves). 



24 



